On the nature of dark energy: the lattice Universe 
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Abstract There is something unknown in the cosmos. 
Something big. Which causes the acceleration of the Uni- 
verse expansion, that is perhaps the most surprising and un- 
expected discovery of the last decades, and thus represents 
one of the most pressing mysteries of the Universe. The cur- 
rent standard ACDM model uses two unknown entities to 
make everything fit: dark energy and dark matter, which to- 
gether would constitute more than 95% of the energy density 
of the Universe. A bit like saying that we have understood 
almost nothing, but without openly admitting it. Here we 
start from the recent theoretical results that come from the 
extension of general relativity to antimatter, through CPT 
symmetry. This theory predicts a mutual gravitational repul- 
sion between matter and antimatter Our basic assumption 
is that the Universe contains equal amounts of matter and 
antimatter, with antimatter possibly located in cosmic voids, 
as discussed in previous works. From this scenario we de- 
velop a simple cosmological model, from whose equations 
we derive the first results. While the existence of the elusive 
dark energy is completely replaced by gravitational repul- 
sion, the presence of dark matter is not excluded, but not 
strictly required, as most of the related phenomena can also 
be ascribed to repulsive-gravity effects. With a matter en- 
ergy density ranging from ~ 5% (baryonic matter alone, 
and as much antimatter) to ~ 25% of the so-called criti- 
cal density, the present age of the Universe varies between 
about 13 and 15 Gyr. The SN la test is successfully passed, 
with residuals comparable with those of the ACDM model 
in the observed redshift range, but with a clear prediction 
for fainter SNe at higher z- Moreover, this model has neither 
horizon nor coincidence problems, and no initial singularity 
is requested. In conclusion, we have replaced all the tough 
problems of the current standard cosmology (including the 
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matter-antimatter asymmetry) with only one question: is 
the gravitational interaction between matter and antimatter 
really repulsive as predicted by the theory and as the obser- 
vation of the Universe seems to suggest? We are awaiting 
experimental responses. 

Keywords Cosmology: theory — Dark energy — Gravita- 
tion — Large-scale structure of Universe 



1 Introduction 

Since the discovery of the cosmic expansion accele ration in 
1998 (e.g. iRiess et aDll998l: Ferlmutter et al.lll999b . one of 
the most debated questions in physics and cosmology has 
been the existence and nature of the so-called dark energy, 
which should account for that unexpected phenomenon. In- 
deed, a repulsive force acting in the Universe space-time de- 
fies any previous physical knowledge, as the only known 
interaction among matter on these large scales is the uni- 
versal Newton-Einstein gravitational attraction. Both the 
classical (Newtonian) and relativistic (Einsteinian) theories 
of gravitation seem to exclude that gravity c an be in some 



way re pulsive. However, in two recent papers IVillatal (120 1 1 



2012ah showed that the general theory of relativity can be 
consistently extended to the existence of antimatter (which 
was unknown at the epoch of the birth of the two theories), 
based on its CPT properties, which imply that matter and an- 
timatter are both gravitationally self-attractive, but mutually 
repulsive. Thus, if our Universe contains a certain amount of 
antimatter (possibly equivalent to the matter one, due to the 
expected matter-antimatter symmetry), the origin of the cos- 
mic speed-up can be easily and naturally explained (together 
with the well-known expansion itself), without any need of 
mysterious ingredients like a physically unknown dark en- 
ergy, or of modifications to the current well-estabUshed the- 
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Thus far, the Universe acceleration, i.e. the presence of 
the elusive dark energy, has been formally ascribed to an 
additional term having a negative pressure in the cosmic- 
expansion equations, in the simplest case corresponding to 
a cosmological constant, perhaps associated to the energy 
of the quantum vacuum. Besides this standard cosmology 
of the ACDM model, various alternatives have been pro- 
posed to explain the accelerated expansion, invoking scalar 
fields or modifications of general relativity, such as exten- 
sions to_extra_dimmsion^^ 

(e.g. lAmendolalboOOl: lOvali et al Jbooot ICarroll et al.ll2004 : 



of this story can be d r awn f r om ma n y works of v arious 



Capozziello et al. 2005 : Napolitano et al. 20121) . For an ex 
tensive and detailed discussion (and bibliography) on alter- 
native models of dark-energ y cosmologies see the recent re- 
view bv lBamba et alj (120121) . 

The main problem with the connection of dark energy 
with the quantum vacuum energy is that the latter would be 
expected to be some 10'^'^ times larger than observed, or 
at least ~ 1 0"^" when consid ering only quantum chromody- 
namics (e.g. Weinberg|l989 ). Recently, starting from the as- 
sun iption of repulsive g ravity between matter and antimat- 



ter, iHaidukovid ( l2012bl) has shown that dark energy could 



be the result of the gravitational polari zation of the quan- 
tum va cuum, while in previous papers (lHaidukovidl201 lb 
2012al) the author showed that also dark matter could be an 
illusion caused by the same phenomenon. This treatment of 
quantum vacuum virtual pairs leads to much more reason- 
able theoretical values for the associated dark energy, but not 
yet in agreement with observational estimates. Another big 
prob lem in this sense is the s o-called "coincidence" problem 
(e.g. [Peebles & Ratrall2003l and references therein): why 
should the dark energy density be so comparable (a small 
factor larger) with the energy density of matter, and why 
just now in the history of the Universe (while, being con- 
stant in the ACDM model, in the past it would have been 
so negligible and in the future it will dominate)? Thus, one 
suspects that actually the repulsive force could be closely 
related to the matter content, as we will show in this paper 
In Sect. 2 we present the basis of general relativity extended 
to the presence of antimatter, while in Sect. 3 we apply it 
in a cosmological model, whose results are submitted to the 
main observational tests in Sect. 4. Conclusions are drawn 
in Sect. 5. 



2 Antimatter gravity 

Antimatter gravity has a long and troubled history, which 
begins with the belief that gravity can be only attractive and 
never repulsive, then this concept passes through phases of 
questioning and contestation, to reach a growing consensus 
towards the possibility that the gravitational interaction be- 
tween matter and antimatter is repulsive. The main steps 



author s e.g. iMorrisonI (Il958l): ISchiffI Jl958[ Il959|): lOood 
(ll96ll ):'Nieto & Goldman' ('l99l'):'Noves & StarsonVl99ll): 



Chardin & Rax (1992): Chardin (1993 1997): Ni (20041): 
Nove j (|2008l)~jCabboletl (l20ld |201 ll): IHaidukovid (I201C 



201 lblaLl2012albl):IVillatal(l201 l[l2012blal) : iBenoit-Levv & Chardinl 
(l2012l) : lDopital(l2012l) : and references therein. 

Here we start from the assumption that antimatter is CPT- 
transformed m atter and, as already partially done in IVillata 



( 1201 iL 1201 2al) . show how gravitational repulsion between 



matter and antimatter is a natural outcome of general relativ- 
ity, and derive the relevant field and cosmological equations. 

2. 1 Equation of motion 

In a metric theory of space-time, the action for a free test 
particle with rest mass m and line element ds is defined as 



S - -m J" di = -m J yjgijvdx'^dx^ . 



(1) 



With x" = t and i'^ - dx''/dt, in terms of the Lagrangian L 
we have 



/ds , 
Ldt, L — —m — =—mJgZJFW. 
dt 



(2) 



With this classical definition of the Lagrangian, the line el- 
ement (or the action) of the test particle, which is a scalar, 
is split into two non-scalar parts, ds/dt (or L) and dt, which 
in particular are both (CP)T-odd, i.e. their CPT-transformed 
counterparts change sign. 

Und er the assump tion that antimatter is CPT-transformed 
matter dViUatalbol Ij), while for matter (M) ds/dt and dt are 
both positive definite, for antimatter (A) they are both nega- 
tive definite. Consequently, we have 



ds 



The canonical momenta are defined as 
BL ^ 



(3) 



(4) 



where m'^ = dx'^/ds. Apart from the sign, they are the three 
spatial components of the covariant version of the energy- 
momentum four-vector. The time component is given by 
the Hamiltonian 

ds 

H = pix' - L = -mgiijU^'x' + m— = mgOf.u'' = -pQ . (5) 

All the terms in Eqs. (|4]i and (|5]i are (C)PT-odd, as expected, 
since four-vectors are always (C)PT-odd, and thus change 
sign from matter to antimatter 
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Through the Euler-Lagrange equation of motion and the 
energy equation, 



dpi _dL AH _ dH 
df ~ dx> ' At ~ dt 



(6) 



which are both (C)PT-invariant with all terms (C)PT-even, 
we can get the well-known geodesic equation, i.e. the four- 
component equation of motion of general relativity. 



d7 ^ " ' 



(7) 



where the Christoffel symbol F^^ represents the (matter- 
generated) gravitation al field. 

As pointed out by Villatal (2011), this equation is com- 
posed of four (C)PT-odd elements. If we CPT- transform all 
the four elements, we obtain an identical equation describ- 
ing the motion of an antimatter test particle in an antimatter- 
generated gravitational field, since all the four changes of 
sign cancel one another. Thus, this CPT symmetry ensures 
the same self-attractive gravitational behavior for both mat- 
ter and antimatter However, if we transform only one of the 
two components, either the field F^^ or the particle (repre- 
sented by the remaining three elements), we get a change 
of sign that converts the original gravitational attraction into 
repulsion, so that matter and antimatter repel each otheiQ. 

2.2 Field equations 

In order to obtain a generally covariant field equation, one 
must derive it from an action composed of a scalar La- 
grangian density £. multiplied by the scalar factor yf^A'^x, 
which guarantee the scalamess of the action itself: 



= /v=5X 



A'^x. 



(8) 



To get the action for a point particle in this form, one usu- 
ally multiplies the time integral in Eq. (|2]i by a space integral 
over the three-dimensional Dirac delta function 5^{x - x(f)), 
where x(f) is the position of the particle at time t. Given a 
system of these particles, its action will be 



^ OT„ 1 6^(x - X„(f)) yjg,,vX^i 



:a^x. 



(9) 



Notwithstanding this somewhat bizarre construction, the 
Lagrangian density maintains its (CP)T-oddness, essentially 
due to As /At - -^/i^^wi^ = !/«"■ And this (CP)T-oddness 



'The geodesic equation for a massless particle, such as a photon, is for- 
mally equal to Eq. 0, except for the parameter i, which can no longer be 
taken as the proper time, being d.v = 0, but it will be an affine parameter 
describing the world line. Thus, a (retarded) photon will be repelled by an 
antimatter-generated gravitational field, and a CPT-transfonned photon, i.e. 
an advanced photon, will be repelled by matter. 



is inherited by the so-called stress-energy tensor through its 
definition in terms of X.- 



dgi''' dx^ dgt'\^ 
which indeed becomes 

6^(x - x„(f)) ««< 



^ S 



^F8 



(10) 



(11) 



The subscript 'S' indicates that we are considering the field 
source, which is (C)PT-odd due to the presence of mJJ, to dis- 
tinguish it from the usual stress-energy tensor T'^'', which is 
(C)PT-even, and that can be obtained from T'^^ by multiply- 
ing by a time integral over the Dirac delta function to get 



y p^(x-X„(.„)) 



(12) 



Thus, while T'"' does not change sign from matter to anti- 
matter, T^^ does, as expected, since an antimatter-generated 
field must be opposite to a matter-generated one. In partic- 
ular, the dominant component for a single non-relativistic 
particle is 



■x(0)^o 



^F8 



(13) 



i.e. essentiall y the gravitational charge of the particle (see 
ViUatall201 11) . which is positive (m'^ = At/As > 0) for mat- 
ter, and negative (m" = At /As < 0) for antimatter On the 
contrary, the time-time component of the usual stress-energy 
tensor, i.e. 



d'^ix-xis)) Q 



r'"^m I -^—=^{u"rAs, 



(14) 



is positive definite, as it must be, representing the energy 
density. 

If the Lagrangian density in Eq. (jSj includes all (matter 
and radiation) contributions to the field source, adding to it 
the (scalar) term for the space-time geometry, one gets the 
so-called Einstein-Hilbert action. 



(15) 



which, through the action principle, yields the Einstein field 
equation 



(16) 



Let us consider an ensemble of matter and antimatter 
point particles where the radiation contribution to the source 
term T^^ is negligible and there is no contribution from a 
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cosmological constant. From Eq. ( fTTT i. we can separate the 
matter and antimatter contributions as 



pressure p measured in the fluid rest frame, and of the fluid 
four- velocity m'^: 



^ S 



SM SA 



Z 



niM- 



^^(X - XM(f)) «M«I 



~ ^ ^ A ■ 



(17) 



Let us now suppose that the number and mass of the mat- 
ter particles are (at least approximately) equal to those of 
the antimatter particles, so as to have equal total masses of 
matter and antimatter If the differences in the d' compo- 
nents are also negligible, what keeps T'^ different from zero 
in Eq. (fTTb (and then prevents the flat space-time of spe- 
cial relativity in Eq. ( fTSl l) is the spatial distribution of the 
various particles, represented by the Dirac delta functions. 
Only if we had all the matter and antimatter particles cou- 
pled two by two in the same space position, we would get 
T'^ - 0, since the gravitational field (or space-time curva- 
ture) produced by each matter particle would be canceled 
by the coupled antimatter particle or, in other words, there 
would be only neutral gravitational charges and no gravity 
at all. But this scenario is highly unlikely, since matter and 
antimatter repel each other (or, alternatively, they would an- 
nihilate each other). 

Starting from a random distribution of particles, one 
can expect that, due to the matter and antimatter self- 
attraction and mutual repulsion, eventually matter and an- 
timatter could be well separate in two distinct regions, pos- 
sibly having first experienced a transient phase where more 
or less massive "islands" of matter are distributed in space 
alternated with similar aggregations of antimatter In this 
phase we would have matter "particles" surrounded by anti- 
matter ones, and vice versa, which means that each particle 
feels more the repulsive effect of the closer opposite parti- 
cles than the attractive one from the more distant like parti- 
cles, which gives rise to a global expansion and could also 
prevent the final aggregation into two distinct blocks. Any- 
way, we focus on this situation, which resembles that of our 
Universe, where matter is organized in supercluster s sepa- 
rated by vast cosmic voids, in which, according to IVillata 
( l2012ah . equivalent amounts of antimatter may be hiddeifl 

In the standard general theory of relativity (i.e. that not 
extended to CPT-transformed matter), the field source in 
Eq. ( fTSl l for an ideal fluid of matter point particles is usu- 
ally expressed in terms of its energy density p and isotropic 



^Previous studies on a matter-antimatter symmetric Universe (e.g. 
ICohen et alj|l99j) seemed to exclude tlie possibility that matter and an- 
timatter domains have sizes smaller than the visible Universe, due to the 
lack of the expected annihilation radiation from the domain boundaries, 
but, unlike in those models, in our scenario annihilation is prevented by 
gravitational repulsion, so that such a lower limit is no longer valid. 



Sfl 



^ fl 



(18) 



This expression can be obtained by taking the average val- 
ues over the particle population of the T'^ components in 
Eq. ( fTTI ). Thus, according to Eq. ( fTTb . we have that Tg^j^ = 
-r^'jj, i.e. also when T^^ is expressed as in Eq. ( fTSl l. As 
a consequence, it is clear that we can not adopt the sim- 
plifying version of the field source in terms of p and p of 
Eq. ( fTSl l when dealing with two different kinds of field- 
generating particles, since in our above model of alternated 
matter and antimatter charges it would be null, totally ignor- 
ing the strong dependence on the particle distribution. We 
can guess that the appropriate source term for such a two- 
charge distribution is proportional to the expression given in 
Eq. (flSl l, which in our model is considered to be equal for 
matter and antimatter (in particular we have pM - Pa - p), 
with rjjjg - = TgJ^g = " Afl ' ^Incc we expect that the 
total effect is repulsive, we set 



^ Sfl 



SAfl 



-a[(p + p)u''u'' -g^^'p].. 



(19) 



just to have a positive. The value of a will depend on the 
specific charge distribution, and we will explore possible 
values in the next section. 

First we check that our guess in Eq. (fl9l l is correct in 
the "generalized-Newtonian" limit (i.e. the Newtonian limit 
extended to the coexistence of positive and negative gravi- 
tational charges). Since in this limit the standard field equa- 
tions ( fT6l ) and ( fTSi ) reduce to the standard Poisson's equation 
V-0 = AnGp, our hypothesized source term in Eq. (fT9) leads 
to a Poisson-like equation of the form 



— -AjiGap , 



(20) 



which we must demonstrate to be appropriate. 

The (generalized) Newtonian potential of a set of point 
particles with mass m„ and gravitational charge m„u^l (mJJ - 
+ 1 for matter and - 1 for antimatter) felt at the space position 



(21) 



Under the assumption that all particles have the same mass 
m„ — m, the Laplacian of Eq. ( f2T1 i at the location x of a given 
matter particle can be expressed by normalizing distances to 
the nearest neighbor one |x - Xi | as 

vV(x) = GamV- — ^ — = -4nGam6\x - xi) , (22) 
|x-xi| 

where a is the dimensionless constant 



" ^ |x- 



x„| 



(23) 
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In a periodic lattice distribution, each elemental cell of size 
~ 2|x-Xi| contains a mass m of matter and a mass m of anti- 
matter, so that md^(x - xi) in Eq. (l22l i represents the density 
p of one of the two components, and Eq. ( l22l i reduces to 
Eq. (Uni, QED. 



3 Cosmological equations 

The electrostatic counterpart of the dimensionless constant 
a is the so-call ed Madelung constant (after the work of 
Madelun3ll918 ), used for the calculation of the binding en- 



ergy of ionic crystal lattice^ Since the infinite summation 
of our "gravitational" a is formally identical to that of the 
Madelung constant, we can refer to the already computed 
values for specific crystal structures. For example, the value 
for the sodium-chloride (rock-salt) structure (i.e. that origi- 
nally calculated by Madelung), which is composed of two 
interpenetrating face-centered cubic lattices, one for each 
of the two ion types, Na^ and CP, is now known to thou- 
sands of decimal digits: lo-NaCil ~ 1.74756. Another simple 
and highly symmetric structure is that of the cesium chlo- 
ride (two interpenetrating cubic lattices resulting in a kind 
of body-centered cubic structure, with each ion at a cube 
center surrounded by eight opposite ions at the vertices of 
the cube). In this latter case the value of the Madelung con- 
stant is lacsCil ~ 1.76267. Two crystal variants of ZnS (zinc 
blende or sphalerite, and wurtzite) give both Iq-zhsI ~ 1.64, 
with a slight difference. When considering also more com- 
plex crystal structures, \a\ is typically found to vary between 
1.5 and 2.5, which we can take as our "trial" range. 

With the Friedmann-Lemaitre-Robertson- Walker metric 
for a spatially homogeneous and isotropic expanding Uni- 
verse with scale factor a{t), the time-time component of the 
field equation ( fTSI l for a matter-antimatter lattice Universe 
(i.e. with the source term given by Eq. ( fT9] l) becomes 



a 4 

- - -7TGa{p + 3p), 
a 3 



(24) 



which clearly shows that in our model, where a was set to be 
positive (as confirmed by Eq. (|23]|). the Universe expansion 
is accelerated (a > 0). In the dust approximation (i.e. no rel- 
ativistic matter or antimatter), the pressure p vanishes, and 
Eq. ( l24l i can be interpreted in generalized-Newtonian terms 
as the (mean) acceleration felt at the surface of a sphere 
of radius a due to the resulting gravitational effect of the 
matter-antimatter content of the sphere. By adding the two 
opposite contributions from matter and antimatter as they 
would be individually (whose sum is equal to zero), Eq. ( l24l i 



This si mile between i onic crystals and Universe structure can also be 
found in lRipaldj tioToh . 



in the dust approximation can be rewritten as 
GMa GMm 



fl = (1 + a)- 



2 ' 



(25) 



where Mm - Mj\ = 4-7Ta^p/3 is the mass/energy of matter 
(and of antimatter) contained in the sphere. Thus, we can see 
that the repulsive contribution is (1 + a) times the attractive 
one. With a in our fiducial range 1 .5-2.5, this ratio is similar 
to that found in the ACDM model between dark energy and 
matter 

A second cosmological equation can be obtained by sub- 
tracting Eq. iTM from the space-space component of the 
field equation (but can also be derived in other ways), to 
eliminate both a and p: 



0' 



-nOap ■ 



(26) 



with < 0, = 0, > for negative, zero, positive spatial 
curvature, respectively. It is evident that in our case A: < 0, 
in contrast with the current standard model for a spatially 
flat Univers^ By rearranging the terms in Eq. (|2FI . one 

gets 



fl2 4 \k\ 
— + -nGap— = — . 
2 3 fl 2 



(27) 



which has the form of the conservation law of the total (ki- 
netic plus potential) mechanical energy per unit mass, equal 
to \k\l2. Using the same trick as in Eq. ( l25T l, we can see that 
also here the ratio between the repulsive potential energy 
and the attractive one is equal to 1 H- o-. From the v = com- 
ponent of the stress-energy conservation law r^.^, - 0, in the 
dust approximation one gets the energy conservation equa- 
tion pa^ - constant. Thus, the potential in Eq. ( |27] | tends 
to zero as fl ^ oo, and a — > Vi^- In this dust-dominated 
phase there is no initial singularity, since a = at a certain 
initial time ti when ai = a(t\) - %nGapt,l?>\k\ (where as usual 
'0' subscripts indicate the present values and we have fixed 
flo = 1). Even though the critical density pcr = 3//~/87rG 
(with H - a I a) has no physical meaning in our model, it 
nevertheless represents a useful and familiar normalization 
for the density parameter Q = plpa- Thus, in terms of //q 
and Qq, Eq. ( [26] ) becomes 



i/5(l +aQo) 1 



■0 



1 -t- oQa 



(28) 



4 Age of the Universe and SN la test 

The ages of the oldest stars in globular clusters con- 
strain the age of the Universe in the range 12-15 Gyr (e.g. 

''Although observations seem to favor <: = 0, a negative ly curved space can 
not be excluded (see e.g. lBenoit-Levv & Chardinll2012h . 
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0.0 0.1 0.2 0.3 0.4 0.5 

Fig. 1 The expansion age of the dust-dominated matter-antimatter Universe as a function of the density parameter CIq for various values 
of a: 1.5 (red), 2 (green), 2.5 (blue); the gray line indicates a possible lower limit of 13 Gyr. 



Krauss & Chabovei 2003). while cosmic microwave back- 



ground (CMB) anisotropy and large-scale structure mea- 



for a flat Universe (e.s. Teamarketal. 2006; 


Jarosik et al. 


2011; 


Komatsu et alj|201 1 ; SuUivan et alj|201 1 


). 



In our dust-dominated, matter-antimatter model, we do 
not find any singularity, so that we can speak of "age of the 
expansion" rather than of "age of the Universe", not know- 
ing when a possible "birth" happened with respect to the 
start of the dust-dominated expansion. However, we may 
suppose that the expansion age is close to what must be com- 
pared with observational results. From Eq. (|28] . the expan- 
sion age, i.e. the time elapsed since the start of the model 
expansion, is 



tfs = to — ti 



1 



Ho VTTffQj 



dfl 



VI - fli/fl 



(29) 



The result of the integration in Eq. (129) is a complicated 
analytic function of aQ.Q. With Hq - 70kms 'Mpc^', in 
Fig. [T]we plot the expansion age as a function of Qq for 
three values of a in the confidence interval: 1.5, 2, 2.5. 
The function starts from l/Z/o ~ 13.97 Gyr at Qo = and 
reaches a maximum of fe,max ~ 15.00 Gyr at Qq ~ 0.092, 
0.069, and 0.055 for a = 1.5, 2, and 2.5, respectively. Then 
it decreases with increasing Qq. As one can see from the 
figure, the expansion age does not yield severe constraints 
on D.O. Indeed, by setting a lower limit to 4 of, e.g., 13 Gyr 
(gray line), the tightest constraint is achieved for a = 2.5 as 



Qo 5 0.25, i.e. nothing particularly unexpected, apart from 
the fact that any lower value (even a baryonic-only Universe 
with Qo ^ 0.05) is not excluded by the oldest-stars con- 
straint. 

Another very important cosmological test is the Hubble 
diagram of Type la supernovae (SNe la), which is so well fit- 
ted by the ACDM model with Qm about 0.27 (Oa ~ 0.73), 
being the historical proof of the expa nsion acceleration . In- 
deed, since the milesto ne works of iRiess et al.l (1 1998b and 
Perlmutter et al.l (119991) . it was evident that distant SNe la 
are fainter than expected from a decelerating Universe. 

In Fig. |2] we plot the correct ed apparent magnit udes of 
the 472 SNe of the data set from lConlev et al.l (1201 11) versus 
their redshift z. For the magnitude correction, ificon - ma + 
a-itis - 1 ) - figC, we adopted the s ame nuisance paramet ers 
found by IConley et all (1201 lb and ISulhvan et all (1201 lb in 
their minimization fits: - 1-37, jSn - 3.18. From 
Eq. (l28T l we derive the luminosity distance in our model as 



1 +z r~ 

VI + aO^) Jo 



At! 



//oVTT^Jo (1 +z')Vl -fli(l +z') 

2(1 +z) r . , r -nO 



//o V 1 + c^^o 



[tanh-' Vl -fli(l 



(30) 



which, with - 70kms 'Mpc^' and a - 2, gives the 
three plotted fits for Qq - 0.05, 0.15, and 0.25. Due to the 
various well-known problem s affecting low-redshift SNe la 
(see e.g. Riess et al.ii2007; .Kessler et al. 20091 : Conlev et al, 




Fig. 3 The same as in Fig. [2] but witii model fits witii a variable as cf oc (red) and a oc a (blue), compared with an o- = constant model 
from Fig. |2] (green) and the ACDM model (dashed). 
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201 ll) . we excluded events with z < 0.02 (36 out of 472, 
cyan symbols in the figure) from the fitting procedure. The 
residuals in the bottom panel show comparably good fits, 
even if Qo = 0.25 yields a better (1.55) due to the 
closer agreement with the large number of low-z SNe. (An 
even x^^^ — 1-53 is obtained for Q.q - 0.28, not 

shown in the figure.) Although we excluded the closest 
SNe, the remaining nearby SNe (say, z < 0.1 ) could still 



be affected by systematic uncertaint ies (see e.g. lKellv et al 



2010HBenoit-Levv & Chardinl2012l) . so that we can not rule 
out low values of Q.q, which actually fit better the most dis- 
tant (z > 1) events, and even a baryonic-only Universe can 
not be excluded. 

On the contrary, if low-z SNe la should not be affected by 
significant bias, it seems that the model fits should be more 
curved (especially at low z) to better match the data. In this 
regard, we recall that our simple model has a main limita- 
tion, since it describes only the dust-dominated phase in the 
Universe history, and this phase would stop in the past when 
a - cii. Even considering the most favorable case Q.q - 0.05 
and ff = 1.5, fli would not be smaller than ~ 0.07, i.e. by far 
too large to meet the epoch producing the CMB radiation. 
We can guess that a must decrease with decreasing time 
while matter becomes hotter, and this would allow the de- 
crease of fli too towards arbitrarily small values. Just for ex- 
ercise and without any serious intention of modeling, we can 
see what would happen with a toy model where a changes 
with time in Eq. ( l28T l as a - aoa^. 

In the Hubble diagram of Fig. |3] we plot the ACDM 
model fit with Q.m = 0.27 compared with three matter- 
antimatter models. The a - constant model (aSlo = 0.3) is 
the same of Fig.|2]with a - 2 and = 0. 15 and is taken as a 
reference for an easier comparison between the bottom pan- 
els of the two figures. We see that with a - a^cfi {[} - 1/2,1 
in our examples) we obtain the above-mentioned curvature 
to better fit the low-z data, similarly to what happens with 
the ACDM model, which indeed fits the data with a ^^^^ 
(1.51) lower than those found for the a = constant models 
(1.53-1.55). Comparable (or better) fits with a variable a 
are obtained with aoQo - 0.67 for a - cq Va C^^^^.^ - 1.51) 
and with ffo^o = 0.85 for a - aoa (jc"^^^ - 1.49). With, 
e.g., Oo = 0.25, in the former case a would vary between 
2.68 and 1.69 from z = to z = 1.5, while in the latter the 
variation range would be a = 3.40-1.36, which appears too 
wide. In any case, as previously found in the a - constant 
models, also here high values of Q.q give lower x^^^'s- How- 
ever, the closeness of these x^^^ values among the various 
models (including the ACDM model) and, even more, the 
existing issue on the low-z data, do not allow us to discrim- 
inate among them on this basis. We notice that all matter- 
antimatter model fits diverge from the ACDM one towards 
high z, predicting significantly fainter SN apparent magni- 
tudes at z > 1 .5. Thus, new data at larger distances, or pos- 



sible revised data for nearby events, should eventually allow 
us to disentangle. 



5 Discussion and conclusions 

In the previous section we have submitted our matter- 
antimatter cosmological model to two specific tests: the 
Universe age and the SN la test, both passed successfully. 
There are other observational constraints that should be 
checked to be in agreement with the model, such as the pri- 
mordial abundances of light elements and the acoustic scale 
of the CMB. Both these constraints have been discussed 



and investigated in detail bv lBenoit-Levv & Chardinl(l2012l) 
in their "Dirac-Milne" cosmology, i.e. a matter-antimatter 
model that appears as a limit case of ours; the one with 
a = 0. Since we have already noticed that in our cosmo- 
logical model a is expected to become very small or null 
in the early stages of the Universe pertinent to primordial 
nucleosynthesis and CMB, we can rely on those results, and 
possibly postpone a detailed study of these issues to future 
works. 

Regarding the radiation-dominated era preceding the 
matter-antimatter dust-dominated one, i.e. when pressure is 
no longer negligible and the energy conservation law be- 
comes pa'^ = constant, one can easily check from Eq. (l26T l 
that also in this case there is no initial singularity, unless 
a ^ with a certain rapidity as f — > fi. Thus, we can 
conclude that in our model no singularity is required, but 
arbitrarily small initial Universe sizes are allowed. Another 
important feature of this model is the absence of the horizon 
problem, since the scale factor acceleration has never been 
negative, and in the earliest stages a can even approach zero. 

Unfortunately, the various tests can not provide strong 
constraints on the value of Qq, but all values between ~ 0.05 
(baryonic-only Universe) and ~ 0.25 (existence of dark mat- 
ter) are possible, even though the higher values seem to be 
favored by lower x^^^'s, in the SN la Hubble diagram fits. In 
any case, in our repulsive-gravity scenario, there seems to be 
no need for mysterious matter in addition to the well-known 
baryonic matter to explain the phenomena for which dark 
matter is usually invoked. Indeed, as shown by lHajdukovic 
( 201 lb, 2012a) . the presence of additional unseen matter at 
galactic scales that would explain, e.g., the galaxy rotational 
curves, can be successfully replaced by the effect of the 
gravitational polarization of the quantum vacuum induced 
by the galaxy mass. On larger, cluster scales, the "observa- 
tion" of potential wells deeper than expected from baryonic 
matter alone, which would allow clustering and would pro- 
duce the we ak gravitational lens ing (see e.g. the recent ob 
servations in 



Dietrich et al ] l2012l) . could just be the effect of 



the presence of surrounding potential "hills" due to antimat- 
ter in the adjacent voids. 
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In summary, starting from the basic assumptions that an- 
timatter is CPT-transformed matter and that our Universe 
is matter-antimatter symmetric, we have developed a cos- 
mological model where, consistently with general relativ- 
ity, gravitational repulsion between the two opposite com- 
ponents is the cause of the accelerated expansion of the Uni- 
verse. This has been done neither with modifications to 
existing well-established theories, nor with the ad-hoc in- 
troduction of unknown entities and dark ingredients. Due 
to the evident absence of matter in the well-known cosmic 
voids, these are the favorite candidates to host antimatter, 
whose invisibility has been discussed and motivated in pre- 
vious works. The resulting lattice structure is well reported 
in the current three-dimensional maps of the observed Uni- 
verse. While in an electrostatic lattice structure (i.e. a crys- 
tal) the alternation of opposite charges (whose interaction is 
attractive) provides a net binding energy in spite of a null 
total charge, the altemance of unlike gravitational charges 
in the cosmos produces a net accelerated expansion in spite 
of the equal amounts of the two components. Similarly to 
the Madelung constant in crystals, the degree of resulting 
repulsive energy is measured by the parameter a, which we 
supposed to be in the range 1.5-2.5, and which multiplies 
the matter (or antimatter) energy density in the cosmological 
equations, thus together providing a single parameter The 
ratio between the repulsive and attractive energies is equal 
to 1 + a, i.e. very close to that found between dark energy 
and matter in the ACDM model, thus solving the coinci- 
dence problem mentioned in the Introduction. In contrast to 
the standard model, the acceleration has never been nega- 
tive and horizon and singularity problems are absent in our 
model. With a in the above confidence interval and with Qq 
in the range 0.05-0.25, the age of the Universe varies be- 
tween about 13 and 15 Gyr. Model fits to the SN la Hubble 
diagram are comparable with that of the ACDM model in 
the observed range z < 1 .4, while they diverge (SNe fainter 
for our model) at higher redshifts, thus offering a future test 
to discriminate between them. Besides dark energy, even 
the existence of dark matter is not needed in our scenario, 
though it is allowed, maybe favored by the SN la test. 

The standard ACDM model is currently the simplest and 
most popular attempt to explain the cosmic acceleration, 
identifying dark energy with the cosmological constant. 
There exists a wide variety of alternative, competing mod- 
els, which are usually more sophisticated but equally com- 
patible with observational constraints. They invoke "dark 
fluids", scalar fields, or geometrical mo difications to genera l 



JOgi 

: alW : 



relativity (see e.g. the recent review bv lBamba et al.ll2012l) . 
However, all of them have the same problem of the simpler 
ACDM model, i.e. no physical justification for the new un- 
known ingredients or geometries, beyond the consistency of 
the models with observational data. 

In conclusion, from the theoretical point of view our 
model appears more elegant and self-consistent than the cur- 



rent dark-energy cosmologies, being based only on well- 
known physical entities and theories, with no need for ad- 
hoc, unknown but dominant, components. Moreover, it 
spontaneously solves several heavy issues like the horizon 
and coincidence problems, the initial singularity, the ap- 
parent matter-antimatter asymmetry. On the other hand, at 
the moment we lack experimental confirmation for the pre- 
dicted repulsive gravity between matter and antimatter, but 
we hope in an answer in a few years f rom the ongoing ex- 
periments (e.g. lKellerbauer et al ] |2008b . 
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